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POPULATION SYNTHESIS OF NEUTRON STARS, STRANGE (QUARK) STARS AND 
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Abstract 

We compute and present the distribution in mass of 
single and binary neutron stars, strange stars, and 
black holes. The calculations were performed using 
a stellar population synthesis code. We follow all 
phases of single and binary evolution, starting from 
a ZAMS binary and ending in the creation of one 
compact object (neutron star, black hole, strange 
star) and a white dwarf, or two compact objects 
(single or binary). We assume that neutron stars 
are formed in the collapse of iron/nickel cores in the 
mass range Mq < M < Mi, quark stars in the range 
Ml < M < M2, and black holes for core masses 
M > M2 and find that the population of quark stars 
can easily be as large as the population of black holes, 
even if there is only a small mass window for their 
formation. 
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1. INTRODUCTION 



Binary population synthesis is a useful tool for study- 
ing the statistical properties of stars, including the 
compact objects. In this paper we wish to address 
the following questions: what is the distribution of 
masses of the compact objects? what is the rela- 
tion between the stellar initial mass and the final 
mass of a compact object, when binary evolution is 
taken into account? given the distribution of com- 
pact object masses, what are the relative numbers of 
different types of objects (neutron stars, quark stars, 
black holes)? 

In section 2 we shortly describe the population syn- 
thesis code used here, in section 2 we present the 
results, and finally in section 4 we present the con- 
clusions. 



2. POPULATION SYNTHESIS CODE 

We use the stellar binary population synthesis code 
consisting of two parts. The single star evolution 
is based on modified formulae from Hurley et al. 
(2000). We have changed the prescription for mass of 
the compact object formed in a supernova explosion. 
We use original Hurley et al. (2000) formulae to ob- 
tain final CO core mass. We use models of Woosley 
(1986) to calculate the final FeNi core mass (for a 
given CO core mass), which will collapse and form 
a compact object during supernova explosion. Fi- 
nally, we include calculations of Fryer and Kalogera 
(2000) to take into account black hole formation both 
through fall back and directly. 

The binary evolution is described in Belczyhski et al. 
(2000). We evolve only binaries where at least one 
star will undergo a supernova explosion. The evolu- 
tion starts at ZAMS. During the course of evolution 
we include the following effects: 

• wind mass loss (standard, Wolf-Rayet, LBV) 

• tidal circularization of binary orbit 

• magnetic breaking 

• conservative/nonconservative mass transfer 

• common envelope evolution 

• rejuvenation 

• hyper-accretion onto compact objects 

• detailed supernova explosion (SN) treatment 

• partial and complete fall back onto compact ob- 
jects in SN 

Many binaries are disrupted in supernova explosions, 
as a result of mass loss and the natal kick. For super- 
nova kicks we use distribution presented by Cordes 
and Chernoff (1997). We continue to evolve all stars, 
both single and binary components, until the forma- 
tion of a stellar remnant. 



2 



Single Stars 



Binary Stars 




Single remnants 



Compact Object 
White Dwarf binaries 



Double Compact 
Object Binaries 



Figure 1. The various routes of compact object formation (compact object = a neutron star, a strange star or 
a black hole). 



3. RESULTS 
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Figure 2. The number per mass interval of neutron 
stars / strange stars / black holes formed along each 
route shown in Fig.l. The solid line corresponds to 
the case of single stellar evolution, the dotted line 
represents the single compact objects formed in bi- 
nary evolution ( group I), the short dashed line shows 
the compact objects in white dwarf binaries (group 
II), and the long dashed line - the compact objects in 
double compact object binaries (group III). 



In Figure 1 we show different formation routes 
of compact objects: neutron stars, strange stars 
and black holes. Compact objects are formed both 
through single and binary stellar evolution. Single 
compact objects are descendants of massive single 
stars but may also be formed as a binary system is 
disrupted in a supernova explosion of one of its com- 
ponents (these are the paths on the left hand side 
of Fig.l). We will denote the single compact objects 
formed as a result of the binary evolution as group 
I. Under favorable conditions some binaries survive 
supernova explosions, and they finally form tight sys- 
tems with compact object/objects. Most of these bi- 
naries will consist of a white dwarf and a compact 
object and the rest will form binaries with two com- 
pact objects (we will denote the compact object in 
binaries with white dwarfs as group II, and the dou- 
ble compact objects as group III). 

In Figure 2 and 3 we show production rates of com- 
pact objects as a function of their final mass for the 
four different formation routes presented in Fig.l. 

Figure 2 shows the number of compact objects per 
mass interval formed along each route of Figure 1. 
We start forming compact objects at mass ~ 1.2Mo 
and their number falls of with the mass of a final 
compact object, as expected for our assumed initial 
mass function ~ M~^-^. A peak around ~ lOAf0 re- 
fiects our relation between ZAMS mass of a progen- 
itor and the final mass of a compact object, which is 
shown for single stars with a solid line in Fig. 4. This 
relation for wide range of progenitor ZAMS masses 
results in a final compact object mass of ~ IOMq. 
This is an effect of a stellar wind which increases with 
a stellar mass, and thus decreases the final mass of a 
compact object flattening out the relation at about 

IOMq. 
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Figure 3. The cumulative fraction of comapct objects (neutron stars, strange stars, black holes) corresponding 
to the differential distribution of Fig. 2. The line types are the same as in Figure 2. 



Figure 3 shows cumulative fraction of compact ob- 
jects as a function of their final mass for different 
formation routes. Note the very fast initial rise of 
cumulative curves for single compact remnants and 
remnants white dwarf binaries. Most of compact ob- 
jects in these groups have mass below 2Mq. Only for 
double compact objects a large fraction of remnants 
have higher final masses. 

In Figure 4 we plotted the mass of a single compact 
object as a function of the ZAMS mass of a progeni- 
tor star. Single compact objects coming from single 
stars are shown with solid line, and the ones com- 
ing from disrupted binary systems are shown with 
points. Points which are over the solid line form a 
group of compact objects whose progenitors gained 
mass in binary interactions, and thus they formed 
more massive compact objects then corresponding 
single stars of the same ZAMS mass. Points which 
arc below the solid lino represent a group of com- 
pact objects whose progenitors were donors in dif- 
ferent mass transfer events, and thus they lost part 
of their mass and formed less massive compact ob- 
jects. There is also a group of compact objects com- 



ing from binary progenitors, which follows exactly 
Mzams-Mfinai relation for compact objects formed 
from single stars. This group comes from wide, non 
interacting binaries, which in the most cases were 
disrupted in a first supernova explosion. The ex- 
ploding component was unaffected by its companion 
so it has left behind a compact object remnant ex- 
actly as it would be formed from a single star. Then 
its companion, also unaffected by binary evolution, 
continues its evolution. Provided that it is massive 
enough, it goes through a supernova explosion leav- 
ing second single compact object remnant, of a mass 
fitting the single star Mzams-Mfinai relation. 



4. CONCLUSIONS 



We have shown the effects of the binary evolution on 
the distribution of masses of compact objects. As ex- 
pected the bulk of the population of compact objects 
have masses below 2M©. While for single stellar 
evolution there exists a unique relation between the 
stellar mass and the compact object mass, there is no 
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such relation when the binary evolution is taken into 
account. The binary evolution works both ways, the 
masses of the compact objects formed in binaries can 
be both smaller and larger than in the corresponding 
case of single stellar evolution. 

The masses of the compact objects are concentrated 
in the low mass range - see Figure 3. Thus even if the 
mass interval (Mi, M2) for the formation of strange 
(quark) stars is small, the fraction of strange stars in 
the whole compact object population can be huge. 
The exact fraction as a function of Mi and M2 can 
be easily read off Figure 3. On the other hand the 
fraction of black holes hardly depends on M2 (the 
cumulative curves flatten above 3 Mq . We conclude 
that the population of quark stars can easily be as 
large as the population of black holes, even if there 
is only a small mass window for their formation. 





Figure 4- The mass of a single compact object as a 
function of the ZAMS mass of the progenitor star. 
The solid line shows the correspondence between the 
initial and final mass for initially single stars. The 
points correspond to compact objects whose progeni- 
tors were initially in binary systems. The three pan- 
els correspond to group I, II, and III, from to bottom 
respectively. Note that binary stars can have more 
then one fate. 
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